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ESTIMATION  OF  NOISE  FIELD  DIRECTIONALITY; 

COMPARISON  WITH  FOURIER  SERIES  METHOD 

INTRODUCTION 

The  possibility  of  estimating  the  directionality  of  a 
stationary  homogeneous  noise  field,  directly  from  the  element 
outputs  of  a  line  array,  was  investigated  in  [1,2,3]  and  found 
feasible  only  for  small  array  sizes,  due  to  ill-conditioning  of 
the  solutions  of  the  fundamental  least-squares  equations  relating 
the  observed  discrete  spatial  correlations  to  the  impingent 
field.  In  addition,  the  possibility  of  incorporating  a  priori 
information  about  the  field  directionality  and  allowing  for 
additive  uncorrelated  noises  at  the  elements  were  considered  in 
[2]. 

Recently,  the  ill-conditioning  associated  with  these 
approaches  was  circumvented  in  [4]  by  employing  a  Fourier  series 
method  for  the  unknown  field.  However,  this  method  requires 
numerous  Bessel  function  evaluations  and  can  have  inaccurate  low- 
order  expansion  coefficients,  leading  to  bias  in  the  estimated 
field  directionality.  Here,  we  will  eliminate  both  of  these 
drawbacks . 
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NOISE  FIELD  CHARACTERIZATION 

Consider  a  (wide-sense)  stationary  homogeneous  noise  field 
characterized  by  directionality  N(f,9,<j>)  that  measures  the  power 
density  spectrum  at  temporal-frequency  f,  due  to  noise  arrivals 
from  direction  (8,<f>).  See  figure  1.  This  characterization 
presumes  that  noise  arrivals  from  different  directions  are 
uncorrelated  and  thereby  precludes  multipath  arrivals,  for 
example . 


Figure  1.  Angular  Geometry 


Suppose  a  collection  of  K  receiving  elements  with  arbitrary, 
but  known,  locations  is  immersed  in  this  field.  The  largest 
amount  of  second-order  information  that  can  be  extracted  from 
these  element  outputs  is  the  set  of  joint  probability  density 
functions  between  elements.  However,  in  so  far  as  estimating 
N(f,9,$)  is  concerned,  there  is  no  need  to  retain  anything  more 
about  the  element  outputs  than  the  (Hermitian)  matrix  G ( f )  of 
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in  order  to  relate  Gkj(f)  to  noise  field  directionality 
N(f,9,+),  consider  the  power  density  spectrum  of  the  elemental 
contribution  due  to  solid  angle  d©  d<|>  sin©  centered  at  (©,$), 


namely 


d©  d$  sin©  N( f , ©, $  )  . 


In  addition,  let  -ck(©,4>)  be  the  time  taken  for  a  noise  arrival 
from  direction  (0,4>)  to  reach  sensor  k  of  a  receiving  array. 
Consequently,  the  transfer  function  applied  to  this  arrival  in 
reaching  the  k-th  sensor  is 


exp[-i2nfrk( ©, +) )  «  Kk(f,©,<fr).  (2) 

Then  the  cross-spectrum  of  the  outputs  of  omnidirectional  sensors 
k  and  j,  owing  to  the  elemental  contribution  (1),  is 

d©  d<fr  sin©  N(f,©,<fr)  Hk(f,e,<M  H_.(f,©,*)*  - 

-  d©  d$  sin©  N{ f , 0, ♦ )  exp[ -i 2 nf { xk ( 9, ♦) (©,♦)}] .  (3) 


Assuming  that  the  noise  arrivals  from  different  directions 
are  uncorrelated,  the  cross-spectrum  Gkj(f)  of  the  outputs  of 
sensors  k  and  j  is  given  by  the  sum  of  components  (3)  over  all 
angular  space: 


G 


ji  n 

J  d©  |  d$  sin©  N  (  f  ,  © ,  $ )  exp[-i2nf  {  Tk(  0,  (  9,  «|>) }  ]  .  (4) 

0  -  n 


This  result  holds  for  1  <  k,j  <  M,  where  M  is  the  number  of 
sensors  in  the  receiving  array.  N(f,0,+)  is  called  the  noise 
fiuld  directionality  at  temporal-frequency  f.  The  product 
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sin©  N(f ,0,4)  could  be  called  the  plane-wave  density. 

Equation  (4),  for  all  k,j,  constitutes  the  totality  of  infor¬ 
mation  about  N( f , 0, 4)  from  the  available  elements.  The  problem 
is  to  estimate  N(f,0,4)  from  measurements  of  {G^j(f)}.  It  should 
be  noted  that  not  all  {G.  .(f)}  contain  independent  information. 
Thus,  Gj^f)  -  G^ff)  for  anc*  G^jff)  ”  G*k(f)  f°r  k»j- 

Also,  for  example,  if  elements  1,2  and  3,4  have  t1-t2  m  T3~T4  ^or 
all  9,4,  (4)  indicates  that  G^tf)  *  <3^  4  (  f )  ?  this  is  due  to  the 
homogeneity  of  the  noise  field. 

Array  processing  techniques,  whether  they  are  standard  delay 
and  add,  weighted,  or  optimum  (adaptive),  result  in  a  preprocess¬ 
ing  of  noise  field  directionality  N(f,0,4).  Attempts  to  then 
estimate  N(f,0,4)  from  these  processed  quantities  must,  in  some 
sense,  undo  what  the  array  processing  has  already  done.  But  why 
should  array  processing  be  used  on  the  elements  at  all,  when  we 
are  interested  in  estimating  the  noise  field  directionality?  An 
"optimum"  estimation  technique  should  accept  matrix  G(f)  as  its 
input  and  emit  an  estimate  of  N(f,0,4)  as  its  output.  This  is 
the  goal  of  this  investigation. 
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LINE  ARRAY 

In  this  report,  we  will  address  only  the  case  where  the  re¬ 
ceiving  array  lies  entirely  on  a  single  line  in  space.  In  figure 
1,  let  9  be  measured  with  respect  to  the  location  of  the  line. 
That  is,  let  the  line  array  lie  along  the  9=0  axis,  namely  the 
z-axis.  See  figure  2.  Then  if  d^  is  the  distance  (measured  down¬ 
ward)  of  the  k-th  element  from  some  reference  point  on  the  line, 
we  have  delay 

dk 

T.  (0,<fr)  =  —  cos9  for  1  <  k  <  M.  (5) 

k  c  - 

The  values  9=0  and  n  correspond  to  endfire  of  the  line  array, 
and  c  is  the  speed  of  propagation. 


t  z 

T 

d, 


Figure  2.  Geometry  of  Line  Array 
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INTEGRATED  DIRECTIONALITY  FUNCTION 


It  follows  from  (4)  and  (5)  that  cross-spectrum 


ji  n 


Gj^f)  »  J  d©  J  d$  sinG  N(f,0,<j>)  exp  -i2nf 


dk  -  di 


cose 


]  -  (6) 


0  -it 


r  dk  '  di 

de  sine  expl-i2nf  - — — cose 


]  n(  f , e) 


for  1  <  k,j  <  M,  (7 


where  we  define  integrated  (or  collapsed  or  averaged)  noise 
di rectionali ty 


N  (  f ,  e ) 


I d* 


N(f ,ef ♦) . 


Notice  that  N(f,0,<fr)  is  defined  in  terms  of  a  coordinate  system 
centered  on  the  line.  According  to  (7),  the  only  quantity  that 
can  be  estimated  about  the  noise  field  directionality  is  the 
integrated  function  N(f,0)  in  (8);  this  is  a  manifestation  of  the 
inherent  conical  symmetry  of  a  line  array  response.  The  problem 
is  to  invert  the  measurements  (G^tf)}  in  (7)  an^  solve  for  the 
quantity  N(f,0),  if  possible. 
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EQUI SPACED  LINE  ARRAY 


For  an  equispaced  line  array,  figure  2  can  be  specialized  to 


d^  =*  k  d  for  1  <  k  <  M, 


(9) 


where  d  is  the  element  spacing.  Then,  (7)  becomes 


Gkj(f)  *  J  d9  sin©  exp^-ia(k-j)  cosej  N(f,9) 


J du  exp £-ia( k-j )uj  N(f,  acos(u) ) , 


(10) 


-1 


where  acos  is  the  principal  value  inverse  cosine  function  and 


_  f  f  c 

“  '  V  o  ■  23' 


U  »  COS0. 


(11) 


The  quantity  fQ  is  the  (design)  frequency  at  which  spacing  d 
would  be  a  half-wavelength: 


2  d. 


(12) 


Also,  u  -  0  corresponds  to  broadside  of  the  line  array. 

At  this  point,  a  change  of  notation  is  very  convenient.  We 
suppress  the  explicit  appearance  of  frequency  f  (it  still  appears 
through  a  in  (11))  and  express  (10)  as  a  spatial  correlation 


C(k-j)  -  J  du  exp( -ia( k- j )u )  B(u), 


(13) 


-1 
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where  we  recognize  that  the  only  dependence  on  k,j  is  through 
their  difference  and  define  the  noise  directionality  function 

f  N(f,  acos(u))  for  -1  <  u  <  1  'i 
B ( u )  -  L  (14) 

l  0  for  | u |  >  1  J 

Finally,  we  modify  spatial  correlation  (13)  to 

1 

C ( p )  -  J  du  exp(-iapu)  B(u).  (15) 

-1 

This  can  be  considered  as  an  integral  equation  for  noise  field 
directionality  B(u),  where  a  is  known,  and  spatial  correlation 
C(p)  -  C  (-p)  is  available  only  for  integer  p  satisfying  |p|  <  M; 
this  was  the  approach  considered  in  (1;  section  2.3]. 
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FOURIER  INTEGRAL  METHOD 


Suppose  that  spatial  correlation  C(p)  in  (15)  were  available 
for  all  continuous  p,  not  just  the  integers  |p|  <  M.  Then  multi¬ 
plying  (15)  by  exp(iopu')  and  integrating  over  all  p,  we  have, 
using  ( 11 )  , 


J  dp  exp(iocpu')  C(p)  -  J  dp  exp(iapu')  J  du  exp(-iapu)  B(u) 


J  du  B ( u )  ~  5 ( u-u ' ) 


B  ( u ' )  if  -1  <  u'  <  1 


0  otherwise 


.  (16) 


That  is, 


B  ( u ) 


J  dp  exp(iaup)  C(p)  for  -1  <  u  <  1. 


(17) 


This  is  an  explicit  integral  relationship  for  the  (integrated) 
noise  field  directionality  B(u)  at  bearing  u,  in  terms  of  spatial 
correlation  C(p)  at  separation  p.  We  will  call  this  the  Fourier 
integral  method  for  the  determination  of  the  noise  field  direc¬ 
tionality. 

As  an  example,  if  the  field  is  composed  of  a  single  plane- 
wave  arrival, 

B ( u )  -  S(u-uQ),  | u0 |  <  1,  (18) 

then  (15)  gives  spatial  correlation 
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C ( p )  -  exp(-iapuQ)  for  all  p,  (19) 

and  (17)  restores  directionality  (18).  We  should  also  note  that 
the  highest  rate  of  variation  of  spatial  correlation  C(p)  is 
exp(+iap),  obtained  when  uQ-+  +1.  If  we  insert  a  higher  variation 
for  C(p)  into  solution  (17),  such  as  (19)  with  | uQ |  >  1,  we  get 
nonzero  values  for  directionality  B(u)  outside  the  allowed  (-1,1') 
range  of  u,  namely  S(u-u  ).  Since  this  is  disallowed  according 
to  (14),  noisy  estimates  of  spatial  correlation  C(p)  require  some 
preprocessing  prior  to  insertion  into  (17);  alternatively,  non¬ 
zero  values  of  B(u)  for  |u|  >  1  might  be  ignored,  but  this  is  not 
an  attractive  approach. 


12 


TR  8599 


INFINITELY-LONG  DISCRETE  ARRAY 


Equation  (17)  gives  the  impression  that  spatial  correlation 
C(p)  is  required  for  all  continuous  p.  But  since  (13)  indicates 
that  C(p)  will  only  be  available  for  integer  p,  we  consider  that 
case  first;  in  fact,  we  consider  spatial  correlation  C(p)  to  be 
known  for  all  integer  p,  which  corresponds  physically  to  an 
inf initely-long  equispaced  line  array.  The  corresponding 
trapezoidal  approximation  to  integral  (17)  for  the  noise  field 
directionality  is  denoted  as 


Bb(u) 


exp(iaun)  C(n) 


for  all  u. 


n  m-to 


But  this  can  be  developed  according  to 


Vu) 


T  w  w 

J  dp  exp(iaup)  C(p)  ^  6(p-n) 


(20) 


n=>-® 


-f-  ou 

u)  ©  ^slu  -  n  ~  -  ^B  U  -  n 


n— • 


n--® 


for  all  u, 


(21) 


where  ©  denotes  convolution,  and  we  used  (11).  Function  Bb(u) 
has  period  2f  /f  in  u;  if  we  define  X  -  c/f,  this  period  is  X/d, 


where  we  used  ( 12 ) . 
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Figure  3.  Approximation  B^(u) 


The  plot  of  approximation  B^u)  to  the  noise  field  direction¬ 
ality  in  figure  3  reveals  aliasing  lobes  separated  by  2fQ/f  on 

the  u-axis.  If  f  <  fQ,  that  is,  if  the  array  is  being  used  below 

its  design  frequency,  then  these  lobes  do  not  overlap,  and  we 
have 

B^(u)  -  B ( u )  for  -1  <  u  <  1.  (22) 

Thus,  exact  recovery  of  the  noise  field  directionality  B(u)  is 
possible  from  knowledge  of  spatial  correlation  C(p)  at  integer  p, 
provided  that 

f  <  fQ,  that  is,  d  <  X/2 ;  (23) 

here,  we  used  X  «  c/f  and  (12).  The  element  spacing  must  be  less 

than  a  half-wavelength  at  the  temporal-frequency  f  of  interest  in 
order  to  avoid  aliasing.  The  discrete  nature  of  the  array  does 
not,  in  itself,  prevent  recovery  of  the  field;  it  is  the  finite 
length  of  a  physical  array  that  causes  problems. 
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If  spatial  correlation  C(p)  contains  a  component  with  too 
high  a  rate  of  variation,  as,  for  example,  (19)  again  with 
| UQ |  >  1,  we  get  3(u)  -  $(u-uo)  as  before.  However,  a  plot  of 
the  corresponding  B^(u)  in  figure  4  (for  uQ  >  1)  reveals  an 
aliased  component  within  the  fundamental  range  (-1,1).  We  would 
be  led  to  believe  that  the  noise  field  directionality  has  a 
component 


6(u  -  u1 ) 


(24) 


which  is  incorrect.  Thus,  the  discrete  nature  of  an  array  can  be 
a  problem  if  the  measured  spatial  correlation  C(p)  contains  dis¬ 
allowed  components,  which  show  up  as  aliased  components  inside 
the  fundamental  range  (-1,1)  of  u.  This  problem  exists  even  if 
the  array  has  infinite  length. 


Figure  4.  Aliasing  of  Disallowed  Component 
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FINITE-LENGTH  DISCRETE  ARRAY 

The  approximate  noise  field  directionality  for  a  finite  array 
of  M  elements  is  a  modification  of  (20): 


BC(U) 


+  ® 


exp( iaun ) 


C ( n )  wn 


for  all  u, 


(25) 


where  weights 

wn=0  for  | n |  >  M.  (26) 

We  consider  that  weights  Wn}  are  samples  of  a  continuous 
function  w(p);  that  is. 


wn  -  w(n) ,  ( 27 ) 

where  function  w  satisfies 

w(p)  -  0  for  | p |  >  M.  (28) 

See  figure  5.  It  then  follows  immediately  from  (25)  and  (20) 
that 
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Bc(u)  -  Bb(u)  ©  W(u),  (29) 

where  window 

+® 

W(u)  -  j|—  J  dp  exp(iaup)  w(p).  (30) 

^  —  OD 

(For  the  case  of  no  weighting,  that  is,  w(p)  1  for  all  p,  then 
W(u)  -  6(u),  and  (29)  reduces  to  (21).)  The  window  W(u)  has 

1  2f 

effective  width  —  —7^  in  u.  (31) 

rl  I 

(For  flat  weighting,  that  is,  w(p)  -  1  for  |p|  <  M,  the  effective 
width  is  half  of  (31);  however,  the  sidelobes  of  W(u)  are  then 
significant. ) 

Although  the  aliasing  of  Bb(u)  in  figure  3  can  be  controlled 
through  satisfaction  of  (23),  the  convolution  result  in  (29) 
reveals  that  the  true  noise  field  directionality  B(u)  will  be 
smeared  by  window  W(u).  This  is  a  result  of  the  finite  length, 
(26),  of  the  array.  Also,  (29)  and  (31)  reveal  that  approxima¬ 
tion  Bc(u)  has  no  super resolution  capabilities;  in  fact,  the 
smaller  that  f  is  chosen  below  design  frequency  fQ,  the  more 
smeared  Bc(u)  becomes.  Thus,  there  is  good  reason  to  operate 
near  the  design  frequency,  that  is,  maximize  f/fQ,  in  order  to 
minimize  the  width  in  (31);  however,  there  is  the  conflicting 
requirement  depicted  in  figure  3  and  (29),  which  points  to 
smaller  values  of  f/fQ.  A  compromise  is  in  order. 
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If  spatial  correlation  C(p)  contains  a  disallowed  component 
such  as  (19)  with  | uQ |  >  1,  (29)  and  figure  4  indicate  the 
presence  of  a  smeared  and  aliased  component  within  the  funda¬ 
mental  range  (-1,1).  Thus,  such  disallowed  components  should  be 
preprocessed  out  of  spatial  correlation  C(p)  before  submission 
into  approximation  (25). 


APPROXIMATE  FIELD 


The  approximate  noise  field  directionality  that  we  will 
consider  at  length,  here,  is  obtained  by  setting  wn  -  1  for 
|n|  <  M  in  (25): 


BJu) 

Cl 


f 

2f 


M— 1 

£ 

n**l-M 


exp(ioun)  C(n) 


for  -1  <  u  <  1 . 


(32) 


This  is  a  sampled  box-car  approximation  to  the  exact  integral 
result  for  B(u)  in  (17).  B_(u)  has  period  2f  _/f  in  u.  If  |C(n)| 

for  | n |  >  M  is  substantially  smaller  than  C(0),  (32)  could  give  a 
good  approximation  to  B(u).  If  not,  then  one  of  the  super¬ 
resolution  techniques,  such  as  maximum  entropy,  could  be  used  to 
effectively  extrapolate  spatial  correlation  C(n)  out  to  n  *  +« 
and  the  transform  carried  out  analytically  in  closed  form. 

Approximations  (32)  and  (25)  have  the  same  form  as  (1;  (41) 
and  (42)],  if  weighting  $(f,u)  there  is  independent  of  u.  Also, 
if  f  -  fQ,  (32)  reduces  identically  to  [1;  (47)  or  (51)]. 
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However,  we  do  not  solve  any  ill-conditioned  equations  here,  when 
f  <  fo;  on  the  other  hand,  we  have  not  minimized  any  error 
measure  either. 

The  actual  numerical  evaluation  of  approximation  Ba(u)  is 
best  done  by  specializing  to  the  particular  argument  values 


B 


a 


M-l 


E 

n-l-M 


exp( i2nkn/K ) 


C  (  n )  - 


■ 

'  M-l 

' 

C(0)  +  2  Re  < 

^  '  exp(-i2nkn/K)  C*(n) 

.  n-1  , 

. 

(33) 


which  can  be  accomplished  by  a  discrete  Fourier  transform.  The 

K  values  of  k,  that  are  swept  out,  cover  an  interval  of  length 

2fQ/f  in  u,  which  is  broader  than  the  length  2  interval  required, 

when  f  <  f  .  See  figure  3. 
o 

The  exp(ix)  function  in  (32)  is  being  sampled  at  increment 

f 

A-au«n-u.  (34) 

x  f 

o 

So,  if  f  <  f  and  |u|  <  1,  then  |4  |  <  n,  meaning  that  exp(ix) 

O  X 

has  at  least  two  samples  per  period.  This  is  well  known  to  be 
the  requirement  for  avoidance  of  overlapped  aliasing  lobes  and  is 
corroborated  by  figure  3.  This  sampling  rate  interpretation  will 
be  very  important  later  when  we  discuss  the  Fourier  series  method 
[4]  . 
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Similarly,  the  spatial  correlation  C(p)  in  (17)  is  being 
sampled  at  increment  A  -  1.  But  since  the  sampling  increment  of 
exp(-ix)  in  (19),  for  basic  elemental  example  (18),  is 


a  u 


f 


(35) 


we  again  have  |  A^  |  <  ji  if  f  <  f  and  |uQ|  <  1.  Thus,  at  least  two 
samples  per  period  are  taken  of  spatial  correlation  C(p),  as 


well,  even  if  values  of  uQ  near  +1  occur. 


DISCUSSION 

It  can  be  seen  from  (34)  and  (35)  that  the  most  troublesome 
cases  will  be  when  frequency  f  is  near  fQ  and  uQ  is  near  +1,  that 
is,  when  the  array  is  employed  at  its  design  frequency  and  when 
arrivals  come  in  near  endfire.  Since  arrivals  come  in  of  their 
own  accord,  no  control  is  had  of  uQ,  except  to  turn  the  line 
array.  And  although  one  could  choose  f  <  fQ  in  order  to 
alleviate  aliasing,  losses  in  resolution  will  then  occur,  as  (31) 
indicates.  Thus,  a  trade-off  is  in  order  in  regards  to  choice  of 
f/fQ;  perhaps,  values  somewhat  less  than  1  are  a  reasonable 
compromise,  as  was  done  in  (1;  figures  3,4,5).  Of  course,  a 
larger  number  of  elements,  M,  always  helps  in  improving 
resolution,  as  shown  by  (31);  this  is  now  a  viable  alternative, 
since  there  is  no  ill-conditioning  as  there  was  in  (1). 
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The  reason  that  we  have  been  able  to  circumvent  the  ill- 
conditioning  is  that  we  have  deferred  the  inherent  sampling  of 
(15),  with  increment  =  1,  until  after  we  solved  the  integral 
equation  for  noise  field  directionality  B(u)  in  explicit  form 
(17).  So,  instead  of  facing  up  to  the  discrete  issue,  as 
explicitly  posed  in  (13),  we  have  put  it  off  as  long  as  possible, 
and  have  then  addressed  it  in  the  various  forms  (20),  (25),  and 
(32),  which  are  reasonable  approximations  to  the  ideal  continuous 
result  (17).  This  procedure  of  temporarily  ignoring  the  discrete 
sampling  associated  with  a  line  array  was  first  presented  in 
[4;  section  3]  in  connection  with  a  cosine  series  expansion  for 
the  field  distribution;  this  latter  procedure  is  fully  developed 
in  the  next  section. 
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FOURIER  SERIES  METHOD 


In  order  to  derive  this  method,  we  return  to  a  combination  of 
(14)  and  (15): 


C(p)  =  J  du  exp(-iapu)  N(f,  acos(u! 


J  d©  exp(-iap  cos©)  sin©  N(f,©)  . 


This  equation  relates  the  measured  spatial  correlation  C(p),  at 
integer  separations  p,  to  the  integrated  noise  field  directional¬ 
ity  N ( f , © )  (see(8)).  We  again  suppress  the  f  dependence  and 
define  plane-wave  density 


A ( 0 )  =  sin©  N ( f , © )  for  0  <  ©  <  n  , 


to  obtain 


C(p)  -  J  d©  exp(-iap  cos©)  A(©) 


As  in  (15)  et  seq.,  spatial  correlation  C(p)  is  known  only  for 
integer  p  satisfying  |p|  <  M.  Function  A(6)  is  the  unknown  field 
function  that  must  be  estimated.  It  is  related  to  the  B  function 
of  the  preceding  section  according  to 


A ( 0 )  -  sin©  B( cos©)  for  0  <  ©  <  n  . 


See  ( 37 )  and  ( 14 ) 
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FOURIER  SERIES  EXPANSION 


We  expand  A(9)  in  a  cosine  series  according  to 


+  08 


A{9)  =  '  a^  cos(q9)  for  0  <  9  <  n 

q-0 


(40 


This  basis  is  a  complete  set  on  interval  (0,n).  See  [5;  page  92]. 
If  we  substitute  (40)  into  (38),  and  interchange  operations,  we 
get 

+«  ji 

C  ( p )  =  y  '  a^  d9  exp(-iap  cos9)  cos(q9)  - 
q*0  ( 


+  a> 


n  aq  ( -i ) q  Jq(otp)  , 

q*0 


(41) 


where  we  used  [6;  9.1.21].  As  a  special  case. 


C(0)  -  n  aQ  ,  (42) 

which  allows  explicit  determination  of  a  . 

o 

Equation  (41)  constitutes  an  infinite  set  of  complex 
simultaneous  linear  equations  for  coefficients  [aq] q  •  If  we 
split  this  equation  into  its  real  and  imaginary  parts,  we  have 
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1  Cr'P> 

+  ® 

+  00 

(-i)q  jq(ap)  -  a2k  (_1)k  J2k(ap)  '  '  (43) 

q-0 

q  even 

k-0 

+  oo 

+  QQ 

I  Ci<P> 

-X>< 

-  T  a2k-l  <-1)k  J2k-l'“P>  • 

q*l 

km  1 

q  odd 

(44) 

We  now  have  two  infinite  sets  of  real  simultaneous  linear 
equations,  one  for  the  even  coefficients,  the  other  for  the  odd 
coefficients,  in  cosine  expansion  (40). 

Since  spatial  correlation  C(p)  in  (38)  is  only  known  for  a 
finite  number  of  discrete  p  values,  namely  integer  |p|  <  M,  there 
is  no  hope  of  solving  for  the  infinite  number  of  unknowns  { }  in 
(43)  and  (44).  What  we  shall  do,  for  the  time  being,  is  to 
ignore  this  limitation  and  pretend  that  C(p)  is  known  for  all 
continuous  p  >  0.  (Of  course,  C(p)  is  then  also  known  for  p  <  0, 
according  to  C(-p)  -  C*(p)  from  (15),  since  field  B(u)  is  real.) 
This  procedure  was  first  propounded  for  the  line  array  in  (4). 

We  multiply  both  sides  of  (43)  by  J2m(ap)/p  and  integrate 
over  p,  to  obtain 


+® 


2m 


( ap) 


+  00 


:r(P)  -  £ 


a2k  '-1’ 


+  CD 


I  dP 


2m 


(  ap ) 


2k 


(  ap) 


k-0 


2m 


(-1) 


m 


2 (  2m) 


for  m  >  1 


(45) 


25 


TR  8599 


here,  we  used  [6;  11.4.6].  In  a  similar  fashion,  from  (44), 


+® 


I  J  dP  ■'-p~P-  Ci'P'  -  E  a2k-l  <-1)k  I  dP  -*=t -  J2k-l(«P> 


+  ® 


+  ® 


J2m-1 * “P1 


k-1 


a2m-l  (-1J 
2 ( 2m-l ) 


m 


for  m  >  1  . 


(46) 


Combining  (42),  (45),  and  (46),  we  have 


ao  -  i  C(0>  ■  i  Cr‘°>  ' 


+  ® 


‘2m  ’  I  '-11"  2">  I  dP  Cr'P» 

0 


a2m-l  *  I  f-1)m  (2”"1>  J  dP  --B=g— ~  Ct 


(P) 


for  m  >  1.  ( 47 ) 


Convergence  of  the  first  integral  at  p  -  0  is  guaranteed  since 
J2m(ap)/p  -*  0  as  p  -*  0,  because  2m  >  2.  The  second  integral  also 
converges  at  p  -  0,  since  J2m-l^aP^  C^(p)/P  *♦  0  as  p  ->  0,  because 
2m-l  >  1  and  C^(0)  «  0.  Thus,  both  integrands  in  (47)  approach 
zero  at  the  origin. 

We  now  have  explicit  integral  relations  for  the  coefficients 
{a^}  in  the  cosine  series  expansion  of  A(9)  in  (40).  They  are 
exact  results  for  {a^},  presuming  that  spatial  correlation  C(p) 
is  available  for  all  continuous  p  >  0.  They  agree  with 
(4;  (13)  —  (15)  ] . 
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EXAMPLE 

Consider  the  same  single  plane-wave  arrival  given  in  (18): 

B(u)  -  S(u  -  uQ)  ,  | uQ |  <  1.  (48) 

Then  (39)  yields 

A( 0)  »  sin©  5(cos0  -  uQ)  -  8(9  -  ©Q)  for  0  <  9  <  it  ,  (49) 

where  0Q  «  acos(uQ),  0  <  0Q  <  it.  Substitution  in  (38)  gives 
spatial  correlation 


C ( p )  «  exp(-iap  cos0Q)  -  exp(-iapuQ)  ,  (50) 


as  in  (19).  When  this  result  is  used  in  (47),  the  coefficients 
are  found  to  be  (see  appendix  A) 


a  ■  —  cos(q  0^)  for  q  >  1 
q  it  o 


(51) 


independent  of  a.  Then  the  summation  on  the  right-hand  side  of 
(40)  becomes  (appendix  A) 


R  ( 0 )  -  [8(9  -  0Q  -  m2  it)  +  5(0  +  0Q 

m—  • 


m2it)  for  all  0. 


(52) 


The  plot  of  this  function  in  figure  6  reveals  that  the  only 
impulse  component  lying  in  the  allowed  range  of  0,  namely  (0,n), 
is  that  at  0  -  0Q.  Therefore, 
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A(  9 )  -  5(0  -  0Q)  for  0  <  0  <  n  ,  (53) 

as  desired.  Thus,  the  use  of  all  the  exact  coefficients  {a^}  in 
(51)  restores  field  distribution  A(0)  precisely. 


Figure  6.  Summation  R(0)  in  (40) 

DISALLOWED  COMPONENT 

The  highest  rate  of  variation  of  spatial  correlation  C(p)  in 

m 

(50)  is  exp(+iap),  just  as  in  (19)  et  seq.  If  we  insert  a  higher 
variation  for  C(p)  into  integral  solution  (47)  for  the  coeffi¬ 
cients,  namely 

C(p)  -  exp(-iapuQ)  ,  uq  >  1  ,  (54) 

we  obtain  [7;  6.693  1&2] 
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independent  of  a.  Use  of  these  coefficients  in  summation  (40) 
gives  reconstructed  field  (right-hand  side  of  (40)) 


(56) 


A  plot  of  this  function  in  figure  7  reveals  that  it  is  spread  out 
over  the  entire  <0,n)  interval;  this  is  in  contrast  to  the 
Fourier  integral  method  in  (19)  which  correctly  restored  a  zero 
field  in  the  fundamental  interval,  namely 


B(u)-5(u-uQ),  uQ  >  1  , 


(57) 
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for  this  example.  Thus,  the  Fourier  series  method  gives  nonzero 
field  values  for  0  <  9  <  re,  even  when  all  the  coefficients  { } 
are  determined  exactly  by  the  integrals  in  (47).  The  use  of 
noisy  estimates  for  C(p)  in  (47)  is  therefore  more  debilitating 
for  the  Fourier  series  method  than  for  the  Fourier  integral 
method,  and  some  preprocessing  (that  is,  low-pass  filtering)  of 
the  available  spatial  correlation  C(p)  values  is  required  prior 
to  insertion  into  (47).  If  this  is  not  done,  a  spurious  back¬ 
ground  will  be  yielded  from  the  Fourier  series  method  in  the 
fundamental  range  0  <  9  <  re,  due  to  "spillover"  from  disallowed 
components  of  C(p). 

Substitution  of  the  reconstructed  field  R(9)  of  (56)  and 
figure  7  into  the  right-hand  side  of  (38)  does  not  restore  the 
spatial  correlation  (54)  for  thic  example  with  uQ  >  1.  JThis  is 
expected,  since  the  coefficients  {a^}  *n  (55)  decay  with  q, 
preventing  summation  (40)  for  A(9)  from  retaining  the  arbitrarily 
narrow  behavior  required  versus  9,  namely  the  delta  function  in 
(49). 
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DISCRETE  ARRAY 

For  an  equispaced  line  array  of  M  elements,  (13)  indicates 
that  spatial  correlation  C(p)  will  be  available  only  for  integer 
|p|  <  M.  We  therefore  adopt,  as  approximations  to  the  exact 
integral  results  in  (47),  the  forms 


*2m  -  I  (-1)m  2” 


M~1  J2m(otn) 

Cr(n) 


2m-l 


\  (-l)m  (2m-l) 


ri-x 

E 

n=*l 

M-l  ,  . 

r — i  J,  1  (  an  ) 

E  —  ci(n 

n-l 


>  for  m  >  1  , 


(58) 


along  with  aQ  -  aQ  -  Cr(0)/n;  see  (4;  (16)— (18) ] .  These  are 
explicit  finite  sums  for  the  approximate  coefficients  { 3  } ,  to  be 
used  in  the  cosine  series  (40)  in  place  of  the  exact,  but 
unknown,  {aq}-  (The  terms  for  n  -  0  in  the  summands  of  (58)  are 
zero  by  virtue  of  the  discussion  under  (47).) 

Several  potential  problems  exist  with  approximations  (58). 
First,  the  increment  in  the  ^(x)  Bessel  functions  in  (58)  is 

f 

A  »  a  -  Jl  —  (59) 

fo 

For  the  design  frequency  f  *  fQ,  (that  is,  d  -  X/2 ) ,  this  incre¬ 
ment  is  n,  which  is  rather  large.  The  plots  of  ^(x)  in  figure 
8,  for  0  <  x  <  20n  and  selected  k  values  between  1  and  35,  reveal 
that  the  low-order  Bessel  functions  are  very  poorly  sampled  at 
values  xn  -  nn,  especially  for  small  n.  For  example,  the  peaks 
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of  J1(x)  at  x  near  .6n  and  1 . 7 u  are  badly  represented  by  samples 
at  x  -  ii  and  2n.  (The  curves  in  figure  8  are  scaled  relative  to 
the  largest  value  of  J^(x)  at  x  near  .6n.) 

Of  more  relevance  are  the  plots  of  weighted  functions  Jk(x)/x 
in  figure  9,  because  this  is  the  actual  integrand  in  (47).  (The 
curves  in  figure  9  are  scaled  individually  for  plotting  appear¬ 
ances,  so  that  each  has  the  same  peak  value.)  We  again  observe 
that  the  low-order  weighted  functions  are  poorly  represented  by 

samples  taken  at  x  *  nit. 

n 

For  larger  arguments  x,  the  "period"  of  J1(x),  indicated 
on  figure  8,  is  approximately  2n;  thus,  we  are  getting  just  2 
samples  per  period  at  &x  -  ji,  which  is  barely  adequate  for  J^(x) 
at  large  argument  values.  For  the  higher-order  weighted  Bessel 
functions,  the  initial  peaks  (near  x  -  k+2)  are  well  represented 
by  samples  at  Ax  -  ji.  In  addition,  the  period  P2  in  figures  8 
and  9  is  greater  than  2n,  for  larger  arguments;  thus,  we  are 
getting  more  than  2  samples  per  period  of  the  higher-order  Bessel 
functions.  Of  course,  eventually,  for  large  enough  x,  all  the 
Jk(x)  have  period  2n.  (See  [6;  9.2.1].) 

These  sampling  considerations  indicate  that  the  low-order 
coefficients  a^  will  likely  be  inaccurate,  especially  when  f  is 
near  fQ,  while  the  higher-order  coefficients  will  not  be  badly 
affected  by  this  particular  feature.  A  numerical  investigation 
of  these  effects  is  undertaken  in  appendix  B.  In  particular, 
coefficients  a^,  a 2,  are  computed  for  a  variety  of  values  of 
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Figure  9.  Weighted  Functions  Jk(x)/x 
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uQ  (cos0Q),  a  ( n f /f o ) ,  M,  and  compared  with  the  exact  values 
a^,  a2,  The  results  quantitatively  confirm  the  above 

expectations . 

The  second  problem  with  approximations  (58)  is  that  the 
increment  in  the  samples  of  spatial  correlation  C(p)  is  -  1. 
The  discussion  surrounding  (35)  is  directly  relevant  again  and 
should  be  reviewed.  Arrivals  near  endfire,  9q  -  0  or  n,  will  be 
most  severely  affected. 

A  third  problem  with  (58)  is  that  M  is  not  infinite;  there¬ 
fore,  the  summands  may  not  have  decayed  sufficiently  to  terminate 
the  summation  at  M-l,  with  negligible  error.  As  seen  earlier, 
for  plane-wave  arrival  (48)-(50),  spatial  correlation  C(p)  does 
not  decay  at  all  with  p,  and  since  [6;  9.2.1] 


vx) 


as  x  -+  +®, 


(60) 


the  integrands  of  (47)  can  decay  very  slowly  with  p. 

Furthermore,  if  a  in  (11)  is  less  than  ji,  J2ni(an)  and 
J2m-l^an^  ^58)  may  not  yet  have  even  reached  their  substantial 
range  of  values  by  the  time  n  reaches  M-l.  To  develop  this 
point,  observe  from  figure  8  that 


J,  (x)  =  0  for  |  x  |  <  k  -  2n. 


(61) 


Therefore 


'J2m^om^  =  0  for  2m  >  an  +  2 n . 
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So,  for  example,  if  M  =  64  and  f  =  f  ,  then  a  -  n,  max  n  -  M-l  - 
63,  and  (62)  indicates  that  the  Bessel  function  is  essentially 
zero  for  2m  >  204.  Thus,  approximate  coefficients  a^,  determined 
from  (58),  will  be  substantially  zero  for  q  >  200;  this  is  the 
limit  that  was  unknown  in  [4;  under  (5b)  and  bottom  of  page 
1651  ]  . 

As  another  example,  if  M  »  64  and  f  *  fQ/2,  then  (62) 
indicates  that  the  approximate  coefficients  a^  for  q  >  105  will 
be  substantially  zero;  this  is  verified  by  [4;  figure  3].  Thus, 
figures  8  and  9,  with  (62),  give  a  quantitative  indication  of 
when  the  Fourier  series  method  will  collapse,  in  terms  of  the 
loss  of  the  higher-order  coefficients  and  the  attendant  degraded 
resolution . 
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RESOLUTION  CAPABILITY 

The  general  situation  is  as  follows:  approximate  coefficient 
a  in  (58)  involves  J  (an)  for  n  *  1  to  M-l.  Reference  to  (61) 

q  q 

therefore  indicates  that  if 


<J  >  »  "max  +  2 “  '  *  <3'.  (63) 

then  a^  =  0.  Thus,  summation  (40)  for  field  A(9)  will  have  zero 
terms  for  q  >  q' .  For  the  same  plane-wave  example  considered  in 
(48)-(53),  this  would  result  in  a  resolution  capability  of  the 
order  of  (appendix  A,  especially  (A-9)  and  (A-10)) 


ji 


n 


(64) 


q ' +  4  a(M-l)  2d( M-l )  2L 

at  broadside,  where  we  used  (63),  (11),  and  defined  L  as  the 
length  of  the  line  array.  However,  the  coefficients  {a^} 
deteriorate  before  q  reaches  q',  typically  for  q  >  q'3/4.  This 
results  in  a  resolution  of  the  order  of 


4  X  2  X 

=  -  -  -  -  .  (65) 

°  3  2L  3  L 


This  is  somewhat  sharper  than  the  standard  quoted  result  of  X/L, 
but  not  significantly  so.  Thus,  the  Fourier  series  method  has 
slightly  better  resolution  than  standard  beamforming,  which 
corroborates  several  of  the  results  in  [4). 
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DISCUSSION 

It  was  demonstrated  in  (51)— (53),  that  for  the  plane-wave 
arrival  of  (48)— (50),  use  of  all  the  exact  coefficients  {a^}  in 
the  Fourier  series  method  restored  the  field  A{9)  precisely  for 
all  0.  However,  when  we  discretize  the  array  and  must  resort  to 
approximations  { }  in  (58),  this  restoration  capability  is  lost, 
even  if  the  array  i s  infinitely  long ;  see  the  tabular  results  in 
appendix  B  for  M  =  100,  1000,  10000,  100000.  This  result  for  the 
Fourier  series  method  is  distinctly  different  from  that  for  the 
Fourier  integral  method,  as  a  review  of  (21)  and  figure  3 
reveals.  In  both  methods,  we  are  presuming  that  f  <  fQ,  that  is, 
that  the  array  is  used  at  or  below  its  design  frequency.  Thus, 
sampling  (in  space)  is  more  detrimental  to  the  Fourier  series 
method  than  to  the  Fourier  integral  method;  this  is  related  to 
the  fact  that  the  latter  employs  a  (single)  Fourier  transform  in 
(17),  whereas  the  former  uses  (numerous)  Bessel  transforms  in 
(47)  . 

The  summation  on  q  in  (40)  for  field  A(0)  cannot  be  carried 
out  to  ®.  However,  when  employed  with  approximate  coefficients 
{a^},  it  should  be  carried  out  at  least  to  the  limit  q'  given  in 
(63),  after  which  {a^}  are  essentially  zero;  this  will  maximally 
preserve  the  resolution  capability  of  the  Fourier  series  method. 
This  procedure  was  not  employed  in  (4;  figures  1,4,6];  thus,  some 
inherent  resolution  of  the  Fourier  series  method  was  lost  in 
those  examples. 
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GRAPHICAL  RESULTS 

We  shall  re-do  the  examples  given  in  [4],  where  the  Fourier 
series  method  was  introduced,  but  now  using  more  coefficients  and 
comparing  the  results  with  the  Fourier  integral  method  presented 
here.  The  first  example  is  that  of  five  plane-waves  with  arrival 
angles  54 0 , 57° , 60 0 , 63° , 66 0 ,  as  given  in  [4;  figures  1  and  2). 
(Angle  90°  corresponds  to  broadside  of  the  line  array. )  The  two 
arrivals  at  57°  and  63°  each  have  twice  the  common  power  of  the 
other  three  arrivals.  The  exact  cosine  series  coefficients  {a^} 
for  0  <  q  <  250  are  plotted  in  figure  10A,  and  are  listed 
numerically  in  table  1  for  0  <  q  <  30.  We  have  normalized  the 
total  power  so  that  the  origin  value  of  the  spatial  correlation 
is  C(0)  *  n;  then  aQ  «  a  ■  1. 

For  a  line  array  with  M  »  64  elements,  employed  at  its  design 

frequency,  f  -  fQ,  the  approximate  coefficients  {a^},  as 

determined  via  (58),  are  given  in  figure*  10B  and  table  1.  A 

comparison  of  the  numerical  results  in  table  1  shows  a  very  large 

discrepancy  between  a^  and  a1  and  between  a2  and  ^2-  However, 

this  discrepancy  decreases  to  about  5%  for  q  -  3  and  4,  and  is 

much  smaller  for  q  >  4.  Comparison  of  the  plots  in  figure  10 

reveals  that  a^  is  substantially  zero  for  q  >  205,  in  agreement 

with  (62)  et  seq. ,  and  that  a  and  a  are  very  similar  for 

q  q 

2  <  q  <  175. 

* 

This  figure  is  not  in  complete  agreement  with  [4;  figure  2]. 
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Table  1.  Cosine  Series  Coefficients 


q 

aq 

§q 

0 

1.0000 

1.0000 

i 

.9977 

.  4837 

2 

-.9906 

-.4407 

3 

-1.9580 

-2.0569 

4 

-.9628 

-.9130 

5 

.9423 

.9265 

6 

1.8349 

1.8428 

7 

.8887 

.8868 

8 

-.8561 

-.8535 

9 

-1.6399 

-1.6395 

10 

-.7806 

-.7784 

11 

.7383 

.7390 

12 

1.3869 

1.3895 

13 

.6463 

.6471 

14 

-.5974 

-.5942 

15 

-1.0938 

-1.0933 

16 

-.4954 

-.4916 

17 

.  4431 

.  4433 

18 

.7809 

.7852 

19 

.3379 

.3376 

20 

-.2857 

'  -.2810 

21 

-.4687 

-.4697 

22 

-.1841 

-.1793 

23 

.1353 

.1335 

24 

.1766 

.1811 

25 

.0433 

.0406 

26 

-.0007 

.0031 

27 

.0790 

.0754 

28 

.0769 

.0794 

29 

-.1114 

-.1157 

30 

-.2857 

-.2851 

For  these  same  parameter  values,  the  reconstructed  field 
distribution,  via  the  Fourier  integral  method  of  (32),  is  depict¬ 
ed  in  figure  11A,  while  that  for  the  Fourier  series  method  of 
(58),  (40),  and  (63),  using  approximate  coefficients  {a^},  is 
displayed  in  figure  11B.  This  latter  figure  is  an  improvement 
over  (4;  figure  1]  for  two  reasons:  250  coefficients  { >  were 
used  instead  of  140,  and  the  angular  sampling  increment  was 
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previously  insufficient  to  track  the  detailed  behavior  of  the 
field  distribution.  The  two  parts  of  figure  11  are  very  similar, 
except  for  the  drift  in  the  Fourier  series  method  near  6  ■  n,  due 
mainly  to  inaccurate  low-order  coefficients  a^  and  a2- 

Both  plots  yield  some  negative  values  for  the  reconstructed 
field  distribution,  due  to  sidelobes  from  the  plane-wave  compo¬ 
nents.  These  can  be  suppressed  at  the  cost  of  decreased  resolu¬ 
tion.  See  (25)— (31)  for  the  Fourier  integral  method.  As  for  the 
Fourier  series  method,  if  approximations  {58)  (to  exact  results 
(47))  used  a  taper,  instead  of  box-car  weighting  out  to  p  -  M-l, 
a  similar  control  of  sidelobes  is  achievable. 

If  all  the  parameter  values  above  are  kept  unchanged,  except 
that  the  arrival  angles  are  squeezed  closer  together,  namely 
56° , 58° , 60° , 62° , 64° ,  we  then  have  the  example  considered  in 
[4;  figures  4  and  5].  The  exact  and  approximate  cosine  series 
coefficients  are  given  in  figures  12A  and  12B,  respectively.  The 
field  distribution  for  the  Fourier  integral  method  is  plotted  in 
figure  13A,  while  that  for  the  Fourier  series  method  is  plotted 
in  figure  13B.  All  five  plane-waves  are  resolved  by  both 
procedures;  in  fact,  the  only  essential  difference  is  the  slight 
drift  in  figure  13B  near  9  -  n,  due  to  poor  values  of  a^  and  §2. 
Figure  13B  is  a  significant  improvement  over  (4;  figure  4],  again 
due  to  additional  coefficients  and  finer  angular  sampling  in  9. 

The  third  example  of  the  Fourier  series  method,  from 
(4;  figure  6],  corresponds  to  three  plane-waves  with  arrivals 
closer  to  endfire,  namely  27°, 30°, 33°.  The  power  level  of  the 
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33°  arrival  is  one-half  of  the  common  power  level  of  the  other 
two  arrivals;  all  other  parameters  are  unchanged.  The  exact  and 
approximate  coefficients  are  given  in  the  two  parts  of  figure  14; 
the  approximation  appears  to  deteriorate  for  q  >  100.  The  field 
distribution  for  the  Fourier  integral  method  is  depicted  in 
figure  15A,  while  that  for  the  Fourier  series  method  is  plotted 
in  figure  15B.  The  major  discrepancy  is  again  the  drift  in  the 
latter  plot  near  0  -  Ft;  this  is  in  spite  of  the  seemingly  poor 
results  for  coefficients  {a^}  in  figure  14B. 

The  final  example  considered  here  is  that  given  in 
[1;  page  15],  namely 


B  ( u ) 


0  for  -1  <  u  <  0  ' 

2u  for  0  <  u  <  1  , 


(66) 


However,  those  earlier  results  were  limited  to  M  <  12  due  to 
ill-conditioning.  The  spatial  correlation  follows  from  (15)  as 


C(p) 


2  2  lexp(-iap)  (1  +  iap) 
a  p  L 


(67) 


The  reconstructed  field  distribution  via  the  Fourier  integral 
method,  for  M  -  64  elements,  is  presented  in  figures  16a  and  16B 
for  flat  weighting  (32),  with  f/fQ  *  1  and  .5,  respectively.  The 
corresponding  plots  for  Hann  weighting,  (25)  and  figure  5,  are 
depicted  in  figure  17.  The  familiar  tradeoff  between  resolution 
and  sidelobes  is  quite  evident.  Perhaps  a  plot  of  both  results, 
with  and  without  weighting,  would  yield  important  information  not 
available  from  either  plot  alone. 


42 


TR  8599 


50  100  15 0 

Figure  10A.  Coefficients  a 
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Figure  10B.  Coefficients  a 


Figure  10.  Coefficients  for  Five  Separated  Arrivals 
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Figure  14B.  Coefficients  a 


Figure  14.  Coefficients  for  Three  Arrivals 
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Figure  16B.  f/fQ  “  *5 


Figure  16.  Directionality  for  Flat  Weighting 
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O  U. 


Figure  17A.  f/fQ  *  1 
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Figure  17B.  f/fo  *  *5 


Figure  17.  Directionality  for  Hann  Weighting 
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SUMMARY 

The  Fourier  integral  method  and  the  Fourier  series  method 
have  very  similar  performance;  the  major  difference  is  the  slow 
drift  in  the  background  level  of  the  noise  field  directionality 
for  the  Fourier  series  method  due  to  inaccurate  low-order  cosine 
series  expansion  coefficients.  There  is  a  rather  large  differ¬ 
ence,  however,  in  terms  of  the  amount  of  computation,  since  the 
Fourier  integral  method  can  employ  a  fast  Fourier  transform  to 
good  advantage,  while  the  Fourier  series  method  requires  numerous 
Bessel  function  evaluations. 

Use  of  the  array  somewhat  below  its  design  -frequency  eases 
the  aliasing  problems  associated  with  both  methods;  but  there  is 
a  tradeoff  connected  with  this  approach,  namely,  a  loss  of  reso¬ 
lution.  Similarly,  weighting  can  be  used  to  suppress  sidelobes, 
but  again,  only  at  the  expense  of  resolution. 

It  has  been  presumed  throughout  this  report  that  the  spatial 
correlation  is  known  exactly,  for  all  required  argument  values, 
without  any  random  error.  In  practice,  the  spatial  correlation 
must  be  estimated  from  a  finite  observation  time  on  random 
processes.  This  limitation  will  further  degrade  the  performance 
of  both  techniques  considered  here;  which  one  will  suffer  most, 
and  by  hew  much,  is  unknown. 
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APPENDIX  A.  EXAMPLE  OF  FOURIER  SERIES  METHOD 


The  single  plane-wave  arrival  was  given  in  (48);  it  yields 
spatial  correlation  (50): 


C(p)  -  exp(-iap  cos90);  0  <  0Q  <  n . 

Substitution  in  (47)  yields  coefficients  [7;  6.693  2&1) 

a2m  "  i  (-1)m  cos(2mS)  1 


2m-l 


\  (-l)ra_1  sin[ ( 2m-l ) 0 ] 


for  m  >  1 


( A-l ) 


( A- 2 ) 


where 

0  -  asin(cos0Q)  -  ^  -  0Q.  (A-3) 

Although  spatial  correlation  C(p)  depends  on  a,  coefficients  {a^} 
do  not.  Also,  observe  that 

cos(2m0)  -  cos ( ran-2m0o )  -  ( — 1 ) m  cos(2m0Q), 

sin[ ( 2m-l ) 0]  -  sin  [  [2m-l]  -  0Q] ]  -  ( -1 )m_1  cos [ ( 2m-l ) 0Q ] ,  (A-4) 


giving 

2 

a„  »  —  cos(q  0_ )  for  q  >  1.  (A-5) 

q  n.  o 

When  these  coefficients,  along  with  aQ  -  1/n,  are  substituted 
in  the  right-hand  side  of  (40),  we  obtain  (letting  0  be 
arbitrary) 
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+® 

R( 9)  a  ^  J  ^  cos ( q@o )  cos(q0)  - 

q-i 


+  ® 


7  +  -  )  [cos{q( 9-0.) )  +  cos(q(0+0  ))]. 
It  II  [_ _ ,  o  o 

q-1 


( A- 6 ) 


But  [8;  page  28] 


+  «D 


+  ® 


+  ® 


2^  +  ^  cos(qt)  -  cos(qt)  =  6(t-m2rt)  ,  (A-7) 

q-1  q*-® 


m»-® 


giving 


+® 

R  ( 0 )  -  [6(0  -  9q  -  m2it)  +  5(0  +  0Q  -  m2n)]  for  all  9.  (A-8) 

m— ® 


This  function  is  discussed  in  (52)  et  seq. 

For  future  reference,  if  the  sum  in  (A-7)  were  terminated  at 
q' ,  we  have 


_1 

2n 


+ 


cos  ( qt ) 


1  sin[ ( 2q '  +  1 ) t/2 ] 
2 n  sin (t/2) 


q-i 

The  first  zero  crossing  of  this  function  is  at 

to  q ' +H  * 

This  is  approximately  the  resolution  of  waveform  (A-9). 


( A-9 ) 


(A-10) 
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APPENDIX  B.  NUMERICAL  INVESTIGATION  OF  (58) 


For  the  single  plane-wave  arrival  given  in  ( 48)— ( 50) ,  the 
spatial  correlation  is 

C(p)  -  exp(-iapuQ)  ,  a  -  n f/fQ  ,  |uQ|  <  1.  (B-l) 

Substitution  in  (58)  gives  approximate  coefficients 

M-l  -  i  v 
r — '  J  •.  (  an ) 

f  *i  -  » -  sin(anu0)  ' 

n-1 


jii 

2 


M-l 

■  n-l 


J  2  (  an  ) 
n  “ 


cos( anuQ) 


n 

2 


( an) 
n 


cosfanu  )  . 
o 


( B-2  ) 


The  exact  coefficients  are  given  by  (51)  as 

£  a„  -  cos ( q0^ )  -  cos(q  acos(u^))  for  q  >  1  ,  (B-3) 

£  CJ  o  o 

and  are  independent  of  a. 

Numerical  values  of  approximations  (B-2)  are  given  in  tables 
B-l,  B-2,  B-3,  respectively,  for  several  values  of  uQ,  a,  and  for 
M  -  100,  1000,  10000,  100000.  The  exact  values,  from  (B-3),  are 
listed  in  the  right-most  column  for  comparison  purposes. 

Several  observations  can  be  made  from  these  tables.  Except 
for  uQ  -  1,  the  sums  in  (B-2)  for  M  -  100  are  not  too  different 
from  what  they  would  have  been  for  M  -  Part  of  this  is  due  to 
the  fact  that  M  -  100  is  considerably  larger  than  the  biggest 
coefficient  order,  10,  that  we  considered  here. 
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The  values  of  the  approximate  coefficients  a^  and  are  poor 
for  a  -  it,  that  is,  for  frequency  f  equal  to  design  frequency  fQ, 
even  for  a  large  number  of  elements  M,  independent  of  arrival 
angle  uQ.  However,  if  a.  is  decreased,  so  that  f  is  well  below 
the  design  frequency  fQ,  d  <  A/2,  then  a^^  and  §2  are  rather  close 
to  a^  and  a2 ,  respectively.  However,  the  loss  in  resolution  is 
unlikely  to  be  tolerable  in  this  case. 

By  contrast,  the  values  of  a^Q  in  table  B-3  are  good 


approximations  to  a^Q,  with 

two 

exceptions : 

uQ  -  1, 

a  m 

n,  all  M  ; 

uo  - 

all 

a,  M  -  100  . 

( B— 4  ) 

That  is,  endfire  arrivals  will  cause  the  most  problems,  as  is 
expected  physically. 
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Table  B-l.  Values  of  ^  a^ 


uo 

a 

M-100 

ih 

M-1000 

f  or : 

M-10000 

M-100000 

ii 

2  al 

m 

all 

0 

0 

0 

0 

0 

.25 

it 

.129232 

.129164 

.129166 

.129166 

.25 

.25 

.  75  it 

.187050 

.187166 

.187169 

.187169 

.25 

.25 

.  5  it 

.223371 

.223390 

.223387 

.223387 

.25 

.25 

.25n 

.244032 

.243504 

.243520 

.243520 

.25 

.5 

ii 

.243202 

.243356 

.243361 

.243361 

.5 

.5 

.  75n 

.371340 

.371278 

.371273 

.371273 

.5 

.5 

.  5  TC 

.446531 

.446295 

.446302 

.446302 

.5 

.5 

.25ri 

.487704 

.487032 

.487014 

.487014 

.5 

.75 

n 

.312112 

.311716 

.311728 

.311728 

.75 

.75 

.  75n 

.548651 

.548352 

.548355 

.548355 

.75 

.75 

.  5  it 

.668571 

.668247 

.668229 

.668230 

.75 

.75 

.  25n 

.730345 

.730503 

.730458 

.730456 

.75 

1 

n 

0 

0 

0 

0 

1 

1 

.75ii 

.675038 

.700277 

.708228 

.710741 

1 

1 

.  5n 

.843507 

.874322 

.884057 

.887135 

1 

1 

.25ii 

.910216 

.953689 

.967453 

.971806 

1 
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Table  B-2.  Values  of  j  §2 


u 

0 

a 

M-100 

f  §2 

M-1000 

for : 

M-10000 

M-100000 

2  a2 

0 

n 

-.774782 

-.774469 

-.774460 

-.774459 

-1 

0 

.  75n 

-.877975 

-.878516 

-.878500 

-.878500 

-1 

0 

.  5n 

-.948405 

-.947527 

-.947499 

-.947498 

-1 

0 

.  25n 

-.984898 

-.987151 

-.987084 

-.987082 

-1 

.25 

n 

-.638731 

-.639064 

-.639054 

-.639054 

-.875 

.25 

.  7  5  n 

-.751196 

-.750903 

-.750906 

-.750907 

-.875 

.25 

.  5n 

-.822242 

-.822027 

-.822058 

-.822057 

-.875 

.25 

.  2  5  it 

-.860117 

-.862071 

-.862055 

-.862057 

-.875 

.5 

n 

-.226752 

-.226437 

-.226428 

-.226427 

-.5 

.5 

.  75n 

-.367826 

-.367496 

-.367515 

-.367515 

-.5 

.5 

.  5  it 

-.444596 

-.445729 

-.445695 

-.445694 

-.5 

.5 

.  25n 

-.485886 

-.486894 

-.486983 

-.486980 

-.5 

.75 

n 

.497761 

.497414 

.497424 

.497425 

.125 

.75 

.  7  5  it 

.275027 

.273883 

.273906 

.273905 

.125 

.75 

.  5n 

.182765 

.181761 

.181708 

.181709 

.125 

.75 

.25ii 

.136903 

.138239 

.138154 

.138149 

.125 

1 

n 

1.992866 

2.080480 

2.108031 

2.116738 

1 

1 

.  75n 

1.105007 

1.155646 

1.171553 

1.176580 

1 

1 

.  5fi 

.969535 

1.031886 

1.051378 

1.057535 

1 

1 

.  25n 

.884033 

.973009 

1.000600 

1.009308 

1 
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Table  B-3.  Values  of  j  a^g 


u 

o 

a 

M-100 

2  S10  for: 

M-1000  M-100 00 

M-100000 

zt 

2  a10 

PU 

IX 

-1.001835 

-1.000047 

-.999998 

-.999996 

-1 

0 

.  75  it 

-.997211 

-1.000083 

-1.000002 

-1.000000 

-1 

0 

.  5n 

-1.004328 

-1.000142 

-1.000004 

-1.000000 

-1 

EM 

.  2  5  it 

-.993657 

-1.000337 

-1.000013 

-1.000003 

-1 

.25 

n 

.818745 

.816853 

.816903 

.816904 

.816895 

.25 

.75/1 

.815431 

.816916 

.816897 

.816895 

.816895 

.25 

.  5n 

.815723 

.817043 

.816890 

.816894 

.816895 

.25 

.25n 

.823767 

.816816 

.816904 

.816892 

.816895 

.5 

n 

-.501787 

-.499985 

-.499935 

-.499934 

-.5 

.5 

.75/1 

-.501841 

-.499903 

-.500002 

-.499999 

-.5 

.5 

.5/i 

-.495030 

-.500171 

-.500005 

-.500000 

-.5 

.5 

.25n 

-.492036 

-.499587 

-.500016 

-.500003 

-.5 

.75 

it 

.589326 

.587351 

.587401 

.587402 

.586426 

.75 

.75/1 

.592357 

.586315 

.586433 

.586429 

.586426 

.75 

.  5/t 

.593192 

.586680 

.586418 

.586426 

.586426 

.75 

.25/1 

.593387 

.586898 

.586445 

.586423 

.586426 

1 

Zl 

1.329489 

1.797572 

1.936303 

1.979867 

1 

1 

.75/1 

.606719 

.882941 

.963235 

.988393 

1 

1 

.  5  it 

.504084 

.856082 

.954934 

.985763 

1 

1 

.25/1 

.252087 

.794347 

.936197 

.979862 

1 
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